Spontaneous Magnetization and Electron Momentum Density in 3D Quantum Dots 
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We discuss an exactly solvable model Hamiltonian for describing the interacting electron gas in 
a quantum dot. Results for a spherical square well confining potential are presented. The ground 
state is found to exhibit striking oscillations in spin polarization with dot radius at a fixed electron 
density. These oscillations are shown to induce characteristic signatures in the momentum density of 
the electron gas, providing a novel route for direct experimental observation of the dot magnetization 
via spectroscopies sensitive to the electron momentum density. 

PACS numbers: 73.22.Dj, 75.75+a, 75.10-b 



Quantum dots, which may be thought of as artificial 
"atoms" , not only provide a nanoscale laboratory for 
fundamental quantum mechanical studies, but also hold 
the promise of novel spintronic and other applications 
[B Q 0- ■ Spin polarization in nanosystems is there- 
fore drawinggreat attention in view of its intrinsic im- 
portance la, |g, [7J , and its possible relationship to other 
phenomena, such as kinks in the positions of Coulomb 
blockade peaks in 2D quantum dots or the so-called 
0.7 conduction anomaly in quantum point contacts 
In 2D dots, experiments point to shell-like filling of lev- 
els in accord with Hund's first rule and a recent spin 
density-functional calculation indicates the presence of a 
ground state with high spin in relatively large dots |l0| . 
In 3D dots, shell structure of energy levels has been re- 
ported and results on the fine structure in the en- 
ergy spectra of gold nanoparticles have been interpreted 
in terms of a ground state spin magnetization larger than 
1/2 However, Ref . fl3l argues that in finite disordered 
systems fluctuations in off-diagonal interaction matrix el- 
ements may suppress ground state magnetization. It is 
clear that spin polarization of quantum dots is a subject 
of considerable current interest and controversy. 

In this work, we develop an exactly solvable many- 
body model Hamiltonian of wide applicability for de- 
scribing correlated electrons confined to a quantum dot. 
The ground state is clearly found to exhibit spontaneous 
magnetization and strong oscillations in spin polariza- 
tion as a function of the dot radius. In order to identify 
observable signatures of the remarkable aforementioned 
oscillations in polarization, we consider the momentum 
density of the electron gas in the dot at some length. The 
key quantity in this connection turns out to be the width 
Ap over which the electron momentum density falls to 
zero around the Fermi cut-off momentum. In particular, 
oscillations in magnetization are shown to induce similar 
oscillations in Ap, indicating that this effect would be 
amenable to direct experimental observation via electron 
momentum density spectroscopies such as angular cor- 
relation of positron annihilation radiation (ACAR) and 



Compton scattering. Our results also bear more broadly 
on the applications of quantum dots in that they point to 
the possibility of obtaining a system in a ferromagnetic 
ground state without requiring low temperatures and/or 
external magnetic fields. 

We consider the Hamiltonian 

The first term (Hq) describes the noninteracting system. 
It corresponds to the spherical square well: V(r) = —Vq, 
for r < R, and V(r) — otherwise, and the associ- 
ated one-particle eigenf unctions tp v with eigenvalues e°; 
v = nlm is a composite index and a labels spin. The 
second term (Hi) describes the interaction and allows 
Coulomb repulsion via the parameter U only between 
electrons of opposite spin — electrons of like spin being 
kept apart by the Pauli exclusion principle. In this con- 
nection, we have carried out extensive numerical esti- 
mates, and find that direct Coulomb matrix elements are 
larger, on the average, by an order of magnitude or more, 
than exchange or other terms for dot radii up to 30 A and 
electron densities ranging from r s — 2 to 5 (Bohr radii) 
[l4j . Spin-orbit effects here are also expected to be small, 
since the confining potential is relatively weak due to the 
screening of the nuclei by the core electrons. For these 
reasons, form (1) should provide a reasonable model for 
the electronic structure of 3D quantum dots. 

The standard one-particle Green function for the many 
electron system can be expanded as Q aa i (rr' , t) = 
J2uw ^ r )¥v( r ')£<7<7'(w',T), with Q aa '(yv',T) = 
— (T T a V(T (T)a v , a ,(0)). It is readily shown that fida^/dr = 
—£,vO,ua — UN-a-a^c, where = e° — fj, (/u, is the chem- 
ical potential), and N- a — Ev' "I' -a fl i/' -u- Then, the 
equation of motion for the Green function Q aa i (yv' , r) is 

d 
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FIG. 1: Typical energy spectrum e„ CT for the interacting 
Hamiltonian (1). Filled (open) circles give majority (minor- 
ity) spin levels, those with equal n joined by solid (dotted) 
lines for clarity. Levels possess a degeneracy of 2Z+1. Horizon- 
tal dashed line marks the highest occupied molecular orbital 
energy, Bhomo- 
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FIG. 2: Spin polarization per particle, £ — (Nf —N{)/N, as a 
function of the dot radius R. Computed values (filled circles) 
are connected by solid lines to guide the eye. Number of elec- 
trons N is indicated on the top axis, (a) Weakly interacting 
case, U — + 0. (b) Strongly interacting case, U = 1.75Utf- 



Since our Hamiltonian possesses the property that 
[Hi, Ho — fiN] — 0, Wick's theorem allows us to 
factorize the preceding ground state average, i.e., 
(T t A_ ct (tK ct (t)4 v ,(0)) = -A_ CT S ff<r ,(iV », where 
N—cr = Ylv' f»' -crt w hh f v i — (j denoting the Fermi occu- 
pation function. A Fourier transform in r immediately 
yields the exact Green function 
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The self-energy thus is S* = UN_ a /h. The inter- 
acting energy levels, t ua , are given by the solutions 
of u> — £ v /h — S* = (frequencies measured with re- 
spect to n/Ti) minimizing the ground state energy, E = 
^^(ej + 5?iS*) f va . Given the number of particles A, 
this defines a set of nonlinear equations for the popula- 
tions, iVf and iVj, of up and down spin states, respec- 
tively. The resulting splitting in energy for states of op- 
posite spin, A = U(N<i — Ni), is uniform, i.e., it does not 
depend on the quantum numbers v. Note that tempera- 
ture enters our formulae only through the Fermi function, 
f va ; all results in the remainder of this article refer to the 
zero temperature limit. 

We comment briefly on the nature of the four essential 
parameters which describe our model: r s , the electron 
density; R, the dot radius or, equivalently, the number 
A of electrons [A = (R/r s ) 3 }; V a , the well-depth; and, U, 
the Coulomb energy. The values of r s and R are related 
to the type and size of the systems considered. In this 
study we use r s — 5, a relatively low density, enabling 
us to consider a wide range of dot radii. Vq may then be 
obtained reasonably in terms of the work function of the 
system studied and the Fermi energy of the free electron 
gas with the same r s . For r s = 5, following Ref. 0, we 
take Vq — 8.62 eV. The choice of the remaining param- 



eter U is somewhat tricky, mainly because correlations 
in quantum dots are not well understood. A handle on 
U can be obtained by considering the average value Utf 
of the Thomas-Fermi (TF) screened Coulomb interac- 
tion, f drdr'\tpv(r)\ 2 VTF(\r — r 'l)l<PM( r ')| 2 j averaged over 
all occupied states in the noninteracting system, with 

/1.56 is the 



vtf(i~) = e exp(—r/lTF)/r, where ^tf = 
TF screening length. However, actual screening in real 
materials is likely to be weaker. For example, in no- 
ble metals, Ref. [l6|, reports an effective screening length 
Ztf/0.73. Using this value, we find an enhanced Coulomb 
interaction given roughly by 1.75 C/tfj which is the value 
we have used throughout this work. 

Fig.njillustrates the nature of the energy levels e vt7 for 
a 12.11 A radius dot, for which N = 96, and U = 27.13 
meV (Utf — 15.50 meV). The self-consistent splitting 
between up and down spin electrons is found to be, 
A = 1.03 eV. The highest occupied molecular orbital 
(HOMO) energy level is, £homo = —5.86 eV, which 
may be thought of as the dot "Fermi energy" . We have 
Af = 67 and Aj, = 29, so that the ground state spin 
polarization per electron is £ = (Af — A^)/A = 0.4. 

We discuss now £ as a function of dot radius with ref- 
erence to Fig. |21 In the weakly interacting limit (U — ► 0), 
A — > 0, still the degeneracy between up and down spin 
in the nl shells is lifted, and as Fig. Eta) clearly shows, 
shells fill in accord with Hund's first rule. Indeed, as the 
number of particles increases with dot radius, the spin 
polarization reaches a peak each time a shell is half filled 
with up-spin electrons, and falls to zero when the shell 
is completed with down-spin electrons. The sequence of 
these so-called "magic numbers", i.e., A values for which 
all occupied shells are completely filled, will of course 
differ here from that in atoms or 2D quantum dots due 
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FIG. 3: Inset: Typical electron momentum density, n(p), and 
the magnitude of its first derivative, n'(p)|; curves are nor- 
malized by V/4iv 3 , where V is the dot volume. The position of 
the peak in |n'(p)| defines the dot "Fermi momentum", while 
its full-width-at-half-maximum (FWHM) defines Ap. Main 
graph: Ap in reduced units of pf (solid line, left scale), and 
spin polarization £ (dotted, right scale) vs. dot radius R. 

to differences in the details of the underlying spectrum. 
The magic numbers up to 106 are indicated on the upper 
axis in Fig.ffia), Turning next to the strongly interacting 
case, Fig. Htb) shows large deviations in £ from a simple 
Hund's rule filling ^t|- This is because A changes with 
each added electron in order to minimize the total en- 
ergy. The first 4 electrons enter the up spin Is and lp 
levels, so that £ = 1, while the next 4 enter the corre- 
sponding down spin levels until C vanishes for N = 8. 
The 9th electron induces a strong change in A, so that 
the Is and lp down spin levels are pushed above the Id 
up spin level, causing the former to empty in favor of the 
latter and the system to be become completely polarized 
again. Similar level reorderings are involved in the other 
"kinks" seen in Fig.|^b). 

In order to identify observable signatures of the re- 
markable changes in the polarization depicted in Fig. [21 
we consider the electron momentum density (EMD), de- 
fined by n(p) = (2-k)- 1 J^duj f(uj)A(p,ui), where / 
is the Fermi function and A is the spectral function 
A(p,w) = -21mG R (p ) uj), with G R {p,u) = G(p,iu n -> 
u> + iS). The inset in Fig. 3 illustrates n(p) and the 
magnitude of its first derivative, |n'(p)|, for the same pa- 
rameter values as in Fig ^ The region of rapid varia- 
tion in n(p) can be characterized via the position, pp, of 
the peak in |n'(p)| and the associated full-width-at-half- 
maximum (FWHM), Ap. In the bulk limit in a metallic 
system pp will tend to the Fermi momentum, where n(p) 
will develop a break and |n'(p)| a (5-function peak, cor- 
respondingly. For simplicity, therefore, we may refer to 
Pf loosely as the dot "Fermi momentum", even th ough 
there are no breaks in the EMD in a finite system (lglll!^ |. 
Fig. |3 considers Ap systematically and shows the pres- 
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FIG. 4: Ap/pF vs. dot radius R, obtained from the ap- 
proximate expression of Eq. (QJ (solid line), compared with 
the exact result (bullets). The baseline corresponding to the 
noninteracting limit, Apo/pF, is given (dotted) for reference. 

ence of a dramatic series of peaks (solid line) spaced more 
or less regularly as a function of R. A comparison with 
corresponding changes in £ (dotted curve), makes it ev- 
ident that peaks in Ap are well correlated with those 
in C, some differences notwithstanding. This correlation 
arises because, at a peak in £, the up and down spin 
EMDs become substantially shifted with respect to each 
other, reflecting the polarization of the system, so that 
the EMD presents a much larger value of Ap overall. In 
short, peaks in Ap clearly are a signature of peaks in f. 

Motivated by the preceding considerations, we have 
obtained an approximate expression for Ap as 

Ap/ PF = Ap /p F + cA(C - Co) R, (4) 

where Co is the spin polarization in the weakly interact- 
ing case (U — > 0), and c is a constant, which depends on 
various dot parameters, with a fitted value in the present 
case of, c = 4.29 x 1CT 2 eV -1 A _1 . The first term in 
Eq. Q describes the U — » limit and can be shown 
to follow the scaling law 0.93 A.(r s /ao)/R The sec- 
ond term incorporates the effect of the spin polarization 
induced through the interaction U . Fig. 0] shows that 
Eq. J3J provides an excellent description of the exact Ap 
data as a function of R. The approximation continues 
to be equally good up to R = 30 A, although for clarity 
results over a shorter i?-range are shown in the figure. 

We emphasize that the form of our solution in Eq. 
does not depend explicitly on that of the noninteracting 
Hamiltonian (Ho). Although details of £ will vary with 
those of the spectrum of Hq, we expect our prediction 
of oscillations in £ with dot radius to be generally ro- 
bust to changes in the shape and dimensionality of the 
confining potential [13, EH- Thus, we expect our results 
to be relevant to almost one dimensional systems. In- 
deed, in metallic nanowires |22]|. Zabala et al. 0] report 
a series of peaks in the magnetic moment per electron 
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as a function of wire radius, correlated strongly with an- 
other observable of the system (the elongation force of 
the nanowire). The physics driving these results is simi- 
lar to that in our work [2^ , despite the fact that Ref . 
considers the stabilized jellium model within the frame- 
work of the spin-dependent density-functional formalism, 
while we treat correlation effects differently, in terms of 
the direct interaction U and neglect various off-diagonal 
exchange and Coulomb contributions 0. 

The relationship between £ and Ap provides a power- 
ful handle for a direct experimental verification of the 
polarization oscillations in quantum dots predicted in 
this study via the use of solid state spectroscopies sensi- 
tive to the EMD, namely, angular-correlation of positron- 
annihilation radiation (ACAR) [positrons have the ad- 
vantage of being local probes, see Weber et al. in Ref.Hif. 
inelastic X-ray scattering (IXS) in the deeply inelastic 
(Compton) regime and scanning tunneling microscopy 
(STM). We have carried out further computations at 
other r s values and find that, generally speaking, po- 
larization effects become less prominent with decreasing 
r s , as the kinetic energy dominates, as well as with in- 
creasing R, as electrons are less confined and U becomes 
weaker. Nevertheless, oscillations in Ap are observable 
even for r s ~ 3, particularly for dot sizes R < 10 A. 

In summary, we have presented an exactly solvable 
model Hamiltonian for discussing the properties of the 
interacting electron gas in the confined geometry of a 
quantum dot. Although the results presented in this ar- 
ticle are based on a confining potential in the shape of a 
spherical square well, computations for other geometries 
and confining potentials would be quite straightforward. 
The ground state is found to exhibit spontaneous mag- 
netization and striking oscillations in spin polarization 
£ as a function of the dot radius R, at the fixed elec- 
tron density considered. The oscillations in £ are shown 
to induce similar oscillations in the width Ap around the 
Fermi cut-off in the momentum density, which we refer to 
for simplicity as the dot "Fermi momentum", providing 
a novel route for direct experimental observation of the 
dot magnetization via spectroscopies sensitive to momen- 
tum density (specially, positron annihilation and Comp- 
ton scattering) . A simple expression for Ap is discussed, 
which gives an excellent approximation to the exact nu- 
merical data as a function of radius. We expect the re- 
sults of the present study to be robust to the details of the 
confining potential and to the approximations inherent in 
the treatment of correlations in our model Hamiltonian. 
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